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1 Introduction 

The slowing down of rotating neutron stars has been observed and modeled for 
decades. The simplest models relate the loss the kinetic rotational energy of the 
star to the emission of magnetic radiation from a rotating dipolar magnetic field 
(MDR), attached to the star jT, |2j [3j [5j. The calculated energy loss by a rotating 

pulsar is assumed proportional to a model dependent power of hi. This relation leads 
to the power law fl = -K fl n where n is called the braking index. The value of n can 
be, in principle, determined from observation of higher-order frequency derivatives 
related to n by [6] 


Jp 

( 2 ) 

When the star is taken as a magnetized sphere, rotating in vacuum , with a constant 
moment of inerta (Mol) and a constant magnetic dipole moment, missaligned at a 
fixed angle to its axis of rotation, n is equal to 3 (for derivation see Section [2]) . 

Extraction of the rotational frequency and its time derivatives from observation 
involves a detailed analysis of the time evolution of the pulses, and of the spectra 
and luminosity of radiation from the related nebulae in a wide range of wavelengths. 
Although data on many pulsars are available in the literature, there are only eight 
pulsars generally accepted to yield reliable data on the pulsar’s spin-down (see Ta¬ 
ble. [H recent compilation [7] and Refs, therein). The third derivative is known only 
for the Crab pulsar [8], and PSR B1509-58 9j. 
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PSR 

Frequency 
(Hz) ' 

n 

Ref. 

B1509-58 

6.633598804 

2.839T0.001 

m 

J1119—6127 

2.4512027814 

2.684T0.002 

HU 

J1846—0258 

3.062118502 

2.65T0.1 

ra 



2.16T0.13 

ra 

B0531+21 (Crab) 

30.22543701 

2.51T0.01 

H3j 

B0540-69 

19.8344965 

2.140T0.009 

mm 

J1833—1034 

16.15935711 

1.8569T0.001 

na 

B0833-45 (Vela) 

11.2 

1.4T0.2 

ra 

J1734—3333 

0.855182765 

0.9T0.2 

mi 


Table 1: Selected pulsars adopted from HUM. 

Examination of Table Q] shows that n = 3 does not agree with observation. There 
have been many attempts to extend/modify the basics of the MDR model. These 
include consideration of magnetic field activity (e.g. [T2J jT8i [19] [20, 21, [22, .23]), su¬ 
perfluidity and superconductivity of the matter within pulsars (e.g. [24; [25] [26]), and 
modifications of the power law and related quantities (e.g. [27U7]). Time dependence 
of the constants in the MDR model has also been considered [18, [28] [29] HQj. In 
particular, time evolution of the inclination angle between spin and magnetic dipole 
axes as been recently addressed 0131!. However, there is no model currently avail¬ 
able which would yield, consistently, the typical spread of values of n as illustrated 
in Table [0 

In this work we focus on determination of the maximum deviation of the braking 
index from the value n = 3 by introducing two modifications of the simple MDR 
model: frequency dependence of Mol, related to the change of shape of a deformable 
star due to rotation, and the macroscopic effect of superfluidity of the pulsar core. 
The correction to the expression the braking index arising from these modifications 
following Glcndenning (32] is derived, and included in the calculation, using four 
realistic Equations of State (EoS) over range of baryonic mass (Mb) . We study 
the relation between the softness of the EoS and the rate of change of the braking 
index as a function of frequency and the M B ). The four EoS’s were also used to obtain 
mass density profiles of the pulsars needed to determine the transition region between 
the crust and core. These results were utilized in simulation of an effect of superfluid 
conditions which eliminate the angular momentum exchange at the threshold between 
the crust and core. The calculation is performed over a full range of frequencies of 
the pulsar from zero to the Kepler frequency and a range of M B from 1.0 to 2.2 M 0 , 
representing the gravitational mass range from about 0.8 to 2.0 M 0 . 
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2 Simple MDR Model 


The total energy loss by a rotating magnetized sphere can be expressed in terms of 
the time derivative of the radiated energy as 

dE 2 2ri 4 • 2 / Q \ 

— = 0 sill a, (3) 

where /! is the magnetic dipole moment of the pulsar, fi — B R 3 . R is the radial 
coordinate of a surface point with the surface magnetic field strength B , D is the 
rotational frequency, and a is the angle of inclination between the dipole moment 
and the axis of rotation [32]. 

Substituting the kinetic energy of a rotating body, dependent on the Mol /, 

e = l -m\ (4) 

into 0 yields 

= -G 2stw “' < 5 > 

Assuming constant Mol, dl/dt = 0, we get 

Cl = — -^-f] 3 sin 2 a. (6) 

Setting K = |^sin 2 a in ([6]) and taking /j, and a constant leads to the commonly 
used braking power law describing the pulsar spin-down due to dipole radiation: 

Cl = -AT] 3 . (7) 

Differentiating (0 with respect to time 

Cl = —3KQ 2 Cl, (8) 


and combining 0 and 0 to eliminate K we get the value of the braking index n 


nil 



= 3. 


(9) 


2.1 MDR model with frequency dependent Mol 

The simple MDR value n = 3 @ is derived taking the /, /i and a as independent 
of frequency and constant in time. However, in reality the Mol of rotating pulsars 
changes with frequency and, consequently with time. [I33([32]. The equlibrium state 
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of a rotating pulsar includes the effect of centrifugal forces, acting against gravity. 
The shape of the pulsar is ellipsoidal with decrease (increase) in radius along the 
equatorial (polar) direction with respect to the rotation axis as pulsar spins down. 
Thus the Mol, and, consequently, the braking index, are both frequency dependent. 
It is convenient to re-write s as 

s( H = - cn4 - < 10 > 

where C = ^/i 2 sin 2 a. Assuming this time that dl/dt is non-zero, differentiation of 
m with respect to time gives 

2 m + ni = -2 cn 3 . (ii) 

Differentiating once more gives 

2m + 2 ni + ni + ni = -6 cn 2 ti. 

Using the chain rule we can write / in terms of D 

dl _ dn dl 
dt dt d£l' 

and obtain 

I = I'(l 
i = ti 2 i" + Tii, 


(13) 

(14) 

(15) 


where the primed notation represents the derivative with respect to 

Substituting the identities shown above into (TIT]) and m, we get the following 
relations for D and D, 


(2i +nr) 

£ _ -6CQQ 2 - D 2 (3/' + QI") 

( 2 / + nr) 

After some algebra it is easy to show that the expression of the braking index as a 
function of angular velocity reads 


(16) 

(17) 


n(D) 


nn _ o ( 3 nr + d 2 /") 

~ (2i + nr) 


(18) 


We note that the magnetic dipole moment of the non-spherical pulsar may in prin¬ 
ciple also change with frequency. Estimation of this effect would require knowledge 
of the origin and distribution of the dipole moment, which is lacking. We therefore 
ignore such change here and restrict ourselves to analysis of the two effects described. 
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3 Calculation Method 


Previous modeling of the braking index using the simple MDR model with constant 
Mol assumed a pulsar with 1.4M© gravitational mass and a radius ~ 10 km. In 
this work, which includes frequency dependent Mol and varying Mb, we solve the 
equations of motion of rotating stars with realistic EoS using two different numerical 
methods. 

3.1 The codes 

The PRNS9 code, developed by Weber |35], (33] , is based on a perturbative approach 
to the equations of motion of slowly rotating near-spherical objects {36, 37J. T° en¬ 
sure the reliability of the PRNS9 code results, we also used the RPN code. This 
code by Rodrigo Negreiros [38] is based on a publically available algorithm, RNS, 
developed by Stergioulas and Friedman [39]. The equations of motion are derived 
directly from Einstein’s equations, following the Cook, Shapiro and Teukolsky ap¬ 
proach [40], described in detail in [41]. Both codes are applicable to rotating stars 
with all frequencies up to the Kepler limit. 

A comparison of the results of the two codes is demonstrated in Figure |T] which 
shows Mol as a function of frequency for a pulsar with the QMC700 EoS and Mb 
= 2.0 M 0 (see Section 1372]) . They differ most, but by less than 10%, as the Kepler 
frequency is approached. The small difference at near zero frequency (about 1.25%), 
due to the difference in behavior of the two low density EoS’s (see Section 13721 ). is 
negligible in the context of calculating neutron star macro-properties. 

3.2 The Equation of State 

An essential input to the calculation of macroscopic properties of rotating neutron 
stars is the EoS. The EoS is constructed for two physically different regimes, the high 
density core and the relatively low density crust. 

The microscopic composition of high density matter in the cores of neutron stars 
is not well understood. We have chosen two EoS’s, which assume that the core is 
made only of nucleons, KDEOvl [42] and NRAPR [43]. These EoS were selected by 
Dutra et al. [44] as being among the very few which satisfied an extensive set of 
experimental and observational constraints on properties of high density matter. In 
addition, we use two more realistic EoS which include in the core the heavy strange 
baryons (hyperons) as well as nucleons. The QMC700 EoS has been derived in the 
framework of the Quark-Meson-Coupling (QMC) model [45] 06} and the Hartree V 
(HV) EoS [47] is based on a relativistic mean-held theory of nuclear forces. The 
maximum mass of a static star, calculated using the Tolman-Oppenheimer-Volkoff 
(TOY) equation, is 1.96, 1.93, 1.98 and 1.98 M 0 for KDEOvl, NRAPR, QMC700 


5 





260 



Figure 1: Mol as a function of frequency for a pulsar with Mb = 2.0 M & as calculated 
with both, RPN and PRNS9 numerical codes. 


and HV, respectively, which is close to the gravitational mass of the heaviest known 
neutron stars [58f ;42j. The EoS’s are illustrated in Figure [2] which shows pressure 
as a function of energy density e in units of nuclear saturation energy density eo = 
140 MeV/frn 3 . We observe that the pressure increases as a function of energy density 
almost monotonically for KDEOvl, NR APR and HV, whereas QMC700 EoS predicts 
a change in the rate of increase at about 4 eo- This change, and the subsequent 
softening of the EoS, happens at the transition energy density marking the threshold 
for appearance of hyperons in the matter. Such a change is not observed in the 
HV EoS. The main reason for the difference between the two hyperonic models is 
that the QMC700 distinguishes between the nucleon-nucleon, and nucleon-hyperon 
interactions (neglecting the poorly known hyperon-hyperon interaction), whereas the 
HV model uses a universal set of parameters for all hadrons. Inclusion of both the 
QMC700 and HV EoS in this work reflects the uncertainty in the theory of dense 
matter in the cores of neutron stars. 


4 Results and discussion 

As detailed in the previous section, the calculation of the frequency dependence of 
the braking index has been done for a multiple combination of codes, EoS’s and Mb 
of the rotating star. We show only typical examples of the results, usually for the 
QMC700 EoS, unless stated otherwise. 
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Figure 2: Pressure vs. energy density e (in units of the energy density of symmetric 
nuclear matter at saturation e 0 ) as predicted by the four EoS’s used in this work. 

4.1 Braking index with frequency dependent Mol 

As a general feature, we find that any appreciable deviation of the braking index 
from the generic value n = 3 is observed only at rotational frequencies higher than 
about 250 Hz. The sensitivity of this deviation to the EoS and Mg is demonstrated 
in Figure [31 As can be seen in Figured the biggest change in the braking index of a 
2.0 M 0 star pulsar is predicted by the HV EoS, followed by the QMC700, reaching ~ 
values 1.75 and 2.15 at 750 Hz, respectively. The two nucleon-only EoS’s, KDEOvl 
and NRAPR behave in a very similar way and predict a larger value of n = 2.5 at this 
frequency. These trends can be directly related to the properties of the EoS’s. Figure 
[3] shows the sensitivity to Mb for the QMC700 EoS. The effect clearly increases with 
decreasing M B . 

4.2 Superfluidity of the core 

The effects demonstrated in Figure [3] were calculated assuming that the whole body of 
a pulsar contributes to the total (core+crust) Mol. However, some theories suggest the 
conditions inside a pulsar are consistent with the presence of superfluid/superconducting 
matter, both in the crust and in the core [241 l25l l26[ [5(7). Superfluid material would 
not contribute to the rotation thus reducing the Mol. 

In this work we considered an extreme case in which the whole contribution of 
the core to the total Mol is removed. This scenario could be realized, for example, if 
either the whole core is superfluid or there is a layer of superfluid material between 
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Figure 3: Braking index as a function of frequency. Left Panel: Braking index cal¬ 
culated for a pulsar with Mb = 2.0 M 0 with all EoS’s adopted in this work. Right 
Panel: Braking index calculated of pulsars with Mb = 1.0 - 2.2 M 0 . 

the core and the inner crust of the star, preventing an angular momentum transfer 
between the core and the crust. Either scenario simply results in removal of the 
contribution of the core to the Mol. 

Elimination of the core contribution can lead to a dramatic lowering of the total 
Mol by more than a factor of three, as shown, as an example, in Figure [4] for the 
1.0 M 0 baryon mass and the QMC700 EoS. The difference between the total and 
crust-only Mol shows a weak frequency dependence with a slight increase above about 
600 Hz. In turn, the reduction of the Mol by removal of the core contribution leads to 
additional changes in the braking index, on top of the changes due to the frequency 
dependent Mol (see Figure [3]) as shown in Figure [U This change is, as expected, 
larger for higher mass stars which contain a more significant proportion of dense core 
material then for lower mass stars being more crust-like throughout. 

4.3 Summary 

The variation of braking index of isolated rotating neutron stars with Mb = 1.0 M 0 , 
1.5 M 0 , 2.0 M 0 and 2.2 M 0 with rotational frequency from zero to the Kepler limit 
within the MDR model with frequency dependent Mol has been investigated. The 
microphysics of the star was included through utilizing realistic EoS’s of the pulsar 
matter. An illustration of the possible effect of superfluidity in the star core has been 
included in the study. 

Compiling results of all models used in this work, including the superfluidity effect, 
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Figure 4: Left panel: Total (crust and core) and the crust-only Mol as a function of 
frequency, calculated for a pulsar with M B = 1.0 M 0 . Right panel: A n represent the 
difference in braking index as a function of frequency between stars with and without 
core contribution to the Mol. Each curve is displayed up to the Kepler frequency of 
the star. 

we deduce a definitive upper and lower limit on the braking index as a function of 
frequency, shown in Figure [5] The maximum change in the braking index is obtained 
with the QMC700 EoS and 1.0 M 0 , the least effect is found for KDEOvl and the 
2.2 M 0 star. Reduction of the braking index from the simple MDR model value 
n = 3 happens only at frequencies that are some significant fraction of the Kepler 
frequency. The calculation predicts that isolated pulsars with the braking index most 
deviating from n — 3 have low M B . For the frequencies of known isolated pulsars with 
accurately measured braking indices (see Table [T]), the reduction away from n — 3 
found in this model, is negligible. 

Finally, we have shown that the simple exclusion of the core due the superfluidity, 
or some superfluid barrier between the crust and core, does not have a strong effect 
on braking in the frequency range of observed isolated pulsars. Further development 
of the idea of a macroscopic description of superfluidity would be interesting. Change 
of the magnetic field due to superfluidity and possible magnetic field expulsion, and 
a consequential increase in surface magnetic field strength B could also be usefully 
explored. 
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Figure 5: Lower and upper limits on values of the braking index as a function of fre¬ 
quency, including results from both numerical codes, all EoS’s, M B , and the superfluid 
condition. The (yellow) shaded area between the two lines defines the location of all 
results within the limits. The pulsar with baryonic M B = 2.2 M 0 and the KDEOvl 
EoS has the highest Kepler frequency and defines the frequency limit in this work. 
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